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PREFACE 


This  Memorandum  was  prepared  for  the  Vela  Analysis  project  under 
the  Advanced  Research  Projects  Agency's  contract  with  RAND.  The  pro¬ 
ject  is  a  broad  study  of  systems  to  detect  above-ground  nuclear  bursts. 
This  Memorandum  discusses  the  possible  use  of  phased  arrays  in  detect1 
ing  the  acoustic  signal  from  atmospheric  nuclear  bursts  and  as  such 
should  contribute  to  ARPA's  R  &  D  program  in  this  field. 


1 


v 


SUMMARY 


An  acoustic  phased  array  serves  to  improve  both  the  signal-to- 
turbulent-noise  ratio  and  the  signal- to-acoustic-noise  ratio.  The 
improvement  in  acoustic  signal-to-noise  ratio  is,  for  isotropic 
noise,  the  array  gain--which  is  calculable  and  is  discussed  herein. 
The  improvement  xn  the  turbulent  signal-to-noise  ratio  cannot  be 
calculated,  in  view  of  the  unknown  correlation  of  noise  in  single 
sensors.  However,  it  is  thought  likely  that  this  improvement  will 
be  at  least  as  much  as  the  acoustic  improvement--especially  at  the 
higher  frequencies  of  greatest  interest  for  small  yields. 

This  Memorandum  considers,  for  simple  unshaded  square  ulane 
arrays,  the  dependence  of  array  gain  upon  frequency,  azimuth,  size 
of  array,  number  of  sensors  in  the  array,  and  uncertainty  in  the 
local  velocity  of  sound.  It  is  shown  that  an  array  of  16  sensors 
at  the  lattice  points  of  a  1/2  x  1/2  mile  square  offers  a  frequency 
dependent  gain  which  rises  to  nearly  10  db  at  0.5  cps  for  any  azi¬ 
muth  and  for  any  sound  speed  in  the  range  of  1000  to  1200  ft/sec. 
Such  an  array  appears  to  be  suitable  for  the  detection  of  acoustic 
signals  from  atmospheric  nuclear  blasts. 
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I.  INTRODUCTION 

The  principal  noise  which  interferes  with  the  reception  of  very 
low  frequency  (below  1  cpsl  acoustic  signals  from  nuclear  bursts  in 
the  atmosphere  ir>  that  which  is  caused  by  local  atmospheric  turbulence 
passing  over  the  microphone.  Such  turbulence,  which  gives  rise  to  a 
fluctuating  pressure ,  leads  to  a  wide-band  noise  output  which  cannot 
be  distinguished,  in  a  single  microphone,  from  noise  which  might  have 
arrived  acoustically.  Even  at  selected  low-wind  sites  this  noise 
from  turbulence  is  many  db  above  true  acoustic  noise,  and  is  many 
score  db  above  theoretical  thermal  noise. 

However,  because  this  turbulent  noise  does  not  propagate  acousti¬ 
cally,  the  turbulent  noises  observed  in  two  microphones  which  are 
separated  by  what  would  be  only  a  fraction  of  an  acoustic  wavelength 
are  usually  little  correlated.  Indeed  the  correlation  is  that  associ¬ 
ated  with  the  scale  of  the  atmospheric  turbulence,  and  much  of  the 
scale  is  less  than  1000  ft,  whereas  t'.e  acoustic  wavelengths  of 
interest  ate  longer  than  1000  f„ .  Consequently,  the  VLF  signal-to- 
turbu lent -noise  ratio  can  be  improved  by  spatial  integration  (i.e., 
by  summing  signals  from  mime  ous  tc- stivers)  . 

i  ite  suRHiing  tab  be  done  in  passive  acoustic  networks  so  as  to 
avoid  the  teat  of  installing  a  large  number  of  microphones  and  their 
adNiH  iatvd  e lei  trim  It h .  In  practice  a  pipe  array  is  used  to  sum 
about  100  signal*  introduced  through  small  holes  distributed  along 
a  line  about  1000  ft  long.  If  the  turbulence  were  wholly  uncor- 
reiated  among  the  100  sample  points,  the  signal-to-turbulent-neise 
ratio  would  be  improved  20  db;  most  of  this  20  db  is  actually  obtained 
over  most  of  the  frequency  range  of  interest. 
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Even  with  this  much  improvement,  the  detection  sensitivity  of 
such  an  acoustic  sensov  is  still  limited  by  nonthermal  noise  many 
score  db  above  thermal  noise.  Consequently,  this  Memorandum  considers 
further  signal-to-noise  ratio  improvements  achievable  by  chased  array 
methods.  For  the  purpose  of  the  following  discussion,  an  individual 
sensor  is  understood  to  consist  of  a  single  microphone  with  its 
associated  pipe  array;  most  of  2,0  db  of  turbulent  noise  suppression 
has  already  been  achieved  in  the  single  sensor,  and  arrays  of  such 

sen.-urs  are  contemplated.  The  individual  sensors  are  assumed  to  be 
omnidirectional  below  1  cps. 
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II.  ARRAY  GAIN 

If  some  number  of  sensors  are  disposed  over  a  region  which  may 
be  several  acoustic  wavelengths  across,  signals  from  these  can  be 
summed  in  order  to  achieve  still  more  improvement  in  the  signal- to- 
turbulent-noise  ratio.  If  signals  from  N  such  sensors  are  summed, 
most  of  10  log  N  db  further  improvement  should  be  achieved,  the  exact 
amount  depending  upon  the  distribution  of  turbulent  scale.  However, 
it  is  necessary  to  introduce  appropriate  time  delays  in  the  individual 
sensor  signals  before  summing,  corresponding  to  the  acoustic  transit 
time  (and  direction)  across  the  array,  and  the  design  of  a  phased 
array  becomes  of  concern. 

Such  an  array  not  only  serves  to  reduce  turbulent  noise,  but  also 
exhibits  gain.  "Gain,"  of  course,  indicates  a  measure  of  the  degree 
to  which  the  phased  array  also  serves  to  discriminate  against  noise 
which  arrives  acoustically  from  distant  sources.  In  practice  the 
signal-to-turbulent-noise  ratio  is  usually  poorer  m  a  single  sensor, 
even  with  a  noise-cancelling  pipe  array,  than  the  signal-to-acoastic- 
noise  ratio--but  probably  not  by  a  great  many  db.  Consequently,  a 
further  improvement  of  the  signal-to-turbulent-noise  rati.o  obtained 
by  a  phased  array  might  be  vitiated  by  acoustic  noise  were  it  not 
that  the  array  gain  also  serves  to  reduce  acoustic  noise.  The  number 
or  db  by  which  a  phased  array  serves  to  impiove  the  signal-to-noise 
ratio  is  not  the  same  for  the  two  noises,  since  one  improvement 
depends  upon  the  distribution  of  turbulent  scale  and  the  other  depends 
upon  the  phasing  properties  of  the  array.  However,  the  two  improvements 
may  turn  out  to  be  comparable,  in  which  case  the  ratio  of  tuibulent 
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noise  to  acoustic  noise  may  be  r.bout  the  same  in  a  phased  array  as 
it  is  in  a  single  sensor. 

In  practice,  the  desired  acoustic  signals  and  the  unwanted 
acoustic  noise  travel  essentially  horizontally  across  the  terrain; 
relatively  little  acoustic  input  arrives  by  steep  angles  from  above. 
Thus  it  would  be  "unfair"  to  count  as  an  improvement  any  array  dis¬ 
crimination  against  noises  arriving  vertically.  For  this  reason 
the  "gain"  discussed  herein  is  normalized  on  radians  rather  than 
on  4tt  steradians--a  departure  from  the  "gain  common!’'  discussed 
in  radar  and  sonar  practice.  Formally,  the  term  "gain"  as  used  here 
is  defined  as  follows; 

Let  there  be,  at  large  distance  at  azimuth  cp,  a  sinusoidal  source, 
frequency  f,  such  that  this  source  produces,  in  a  single  sensor,  an 
output  of  unit  amplitude.  Let  S(0,  cp,  c,  c  ,  f)  be  the  output 
amplitude  from  the  whole  phased  array  of  N  sensors,  caused  by  this 
same  source,  when  the  array  is  phased  to  azimuth  8  on  the  assumption 
that  the  local  sound  speed  is  c^  although  the  correct  value  is  c. 

Then 

S2(9,  9,  c,  c  f) 
gain=  - - - 

2tt  J  s2(0’  <P»  c>  V  f>dcP 

o 

“  a  function  of  N,  the  array  geometry, 

9,  f,  c,  and  c 

o 


Ill .  ARRAY  CONFIGURATION 


An  array  might,  of  course,  be  assembled  in  any  configuration, 
and  perhaps  certain  particular  arrangements  would  be  advantageous. 
Indeed,  if  all  azimuths  are  equally  interesting,  then  circular  arrays 
would  seem  to  be  indicated.  However,  real  estate  isn't  ordinarily 
available  in  circular  patches.  Consequently,  arrays  disposed  in  a 
square  region  have  been  chosen  for  examination. 

For  the  purpose  of  illustration  arrays  disposed  in  a  1/2  x  1/2 
statute  mile  square  are  considered.  This  seems  to  be  a  reasonable 
size  for  detecting  nuclear  bursts.  However,  it  should  be  noted  that 
the  curves  shown  below  can  he  shifted  to  any  other  size  of  square 
merely  by  shifting  the  frequency  scale  so  as  to  preserve  the  geometry 
in  wavelength  measure. 

Of  the  many  possible  arrangements  of  sensors  within  this  square. 

consideration  is  limited  to  the  set  of  4,  9,  16,  25 . sensors 

disposed  on  the  lattice  points.  This  set  may  not  be  opt: mu. 
especially  if  it  should  be  desirable  to  discriminate  specifically  in 
favor  of  or  against  certain  fixed  azimuths,  but  it  serves  adequately 
to  illustrate  the  amount  of  array  gain  achievable  with  various 
numbers  of  sensors.  Inasmuch  as  gain  is  of  interest,  rather  than 
side  lobe  control,  shading  schemes  are  not  considered. 
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IV.  GAIN  VERSUS  FREQUENCY  FOR  SEVERAL  ARRAYS 
The  limiting  case  of  this  set  of  arrays  would  be  a  continuously 
sensitive  square,  consisting  of  an  infinite  number  of  infinitesimally 
small  sensors,  each  one  properly  phased  and  then  all  summed.  This 
case  furnishes  the  limit  which  various  arrays  of  finite  numbers  of 
sensors  approach,  and  it  serves  as  a  useful  reference  against  which 
to  intercompare  the  other  arrays. 

The  solid  curve  of  Fig.  1  shows  the  array  gain  (in  db)  versus 
frequency  for  such  a  continuously  sensitive  1/2  x  1/2  mile  square 
phased  to  0  azimuth  (throughout  the  Memorandum  azimuths  are  measured 
from  the  normal  to  one  edge  of  the  square).  Of  course,  the  directivity 
pattern  of  the  square,  and  thus  the  array  gain,  varies  with  the 
azimuth  to  which  it  is  phased;  the  dotted  curve  in  Fig.  1  shows  the 
gain  versus  frequency  for  the  same  square  array  phased  to  45°.  It 
is  seen  that  the  two  curves  differ  negligibly  (the  variation  of  gain 
with  azimuth  is  greater  for  the  finite  arrays  considered  below).  In 
Fig.  1  it  is  assumed  that  the  velocity  of  sound  is  uniform  over  the 
array  and  is  known  exactly  (c  -  Cq  =  1100  ft/sec),  so  no  error  in 
phasing  occurs.  The  solid  curve  of  Fig.  1  represents  a  more  or  less 
Ideal  case;  it  is  repeated,  for  reference,  in  Figs.  2,  2,  and  4. 

In  practice  the  local  velocity  of  sound,  c,  is  not  likely  to  be 
known  exactly,  so  as  to  permit  exact  phasing.  Indeed,  it  would  be 
expensive  and  difficult  to  measure  it  and  to  adjust  the  phasing  net¬ 
work  as  c  changes.  It  would  be  preferable  to  be  able  to  assume  a 
nominal  local  velocity,  Cq  ,  phase  the  signals  on  this  basis,  and  then 
accept  the  degradation  oc  gain  which  may  result  if  the  true  velocity 


nuously  sensitive  array 
1/2  mile 


Frequency  (cps) 


of  suund  is  different  from  nominal.  Consequently,  it  is  not  only 
the  ideal  gain,  when  there  is  no  error  in  sound  speed,  which  is  of 

interest  here,  but  also  the  sensitivity  of  gain  to  error  in  sound 
speed . 

Furthermore,  the  dependence  of  gain  upon  frequency  may  be  quite 
different  at  different  azimuths,  and  so  too  may  the  dependence  of 
gain  on  error  in  sound  speed.  There  would  be  little  merit  in 
achieving  a  high  gain  at  some  particular  azimuth  if  the  same  array 
yielded  much  poorer  gain  at  another  azimuth  or  at  a  slightly  differ¬ 
ent  sound  speed.  Consequently,  the  envelope  of  maximum  and  minimum 
gain  achieved,  versus  frequency,  over  the  range  of  all  possible 
azimuths  and,  simultaneously,  over  a  range  of  sound  speeds  is  of 
greatest  interest.  It  is  this  envelope  which  is  shown  in  Figs.  2, 

3,  and  4  for  arrays  of  4,  9,  and  16  sensors.  On  each  of  these  the 
solid  curve  of  Fig.  1  is  repeated,  and  Figs.  3  and  4  also  repeat, 
dotted,  the  envelopes  shown  in  the  figures  which  proceed  them;  thus 

Fig.  4  permits  an  intercorapar ison  of  all  the  curves  from  the  previous 
figures . 

To  construct  these  curves,  the  arrays  were  assumed  to  be  phased 
on  the  basis  of  a  local  sound  speed  cq  -  1100  ft/sec.  The  gain  was 
computed,  at  each  frequency,  for  every  combination  of  c  «  1000(20) 
1200  ft/sec,  azimuth  =  0(1)45  degrees;  at  each  frequency  the  maximum 
and  minimum  gains  are  plotted.  The  mathematical  expressions  involved 
are  given  in  the  Appendix. 
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— DISCUSSION  OF  ARRAY  GATNg 

For  c  »  c  ,  and  at  any  one  azimuth 

o  y  azimuth,  the  gams  of  the  discrete 

arrays  begin  at  0  db  at  low  frequency  curv, 

quency,  curve  upward  to  a  slope  of 

approximately  3  db/octave,  rare  or  less  following  the  t,™  for  the 

continuously  sensitive  array,  and  then  roll  over  to  approach  their 
respective  high-frequency  asymptotic  values  (gain  .  „) .  It  „ 

interesting  that,  in  the  low-  and  middle-frequency  regions,  the  fewer 
the  number  of  sensors,  the  higher  the  gain. 

These  curves  approach  their  asymptotes  quite  slowly,  even  at  a 

single  azimuth,  and  with  no  error  in  sound  speed,  the  gain  fluctuates 

quite  erratically  about  the  asymptote.  The  expressions  for  gain 

contain  Bessel  functions  of  order  zero-  t-h1ic  u, 

oraer  zero,  thus  these  terms  tend  to  die 

out  asymptotically  as  (frequency)"172. 

For  these  discrete  arrays  it  is  seen  that  the  variation  of  gain 

with  azimuth  and  with  sound  speed  is  slight  up  to  a  frequency  at 

which  the  envelope  abruptly  blooms  open  to  a  rather  extreme  ■ 

lamer  extreme  variation. 

The  lower  boundary  of  gain  reaches  „  m 

g  reaches  a  maximum  when  the  distance  between 

sensoi  about  0.4  wavelengths  (about  0.33  >.  for  4  sensors),  and  at 
higher  frequencies  this  lower  boundary  varies  erratically  downward. 
(Because  of  the  error  In  phasing  caused  by  error  in  sound  speed,  the 
lower  boundary  goes  to  negative  gain  at  higher  frequencies.)  If  the 

frequency  at  which  the  lower  boundary  is  a  maximum,  is  adopted  as 
the  maximum  usable  frequency,  then: 


w  ...  Minimum  Gain 

l  Maximum  Usable  Freouencv  »•„ 

Number  of  Scnsnrq  ,  .  quenc>  at  Maximum 

'  °  *  Usable  Frequency  (db 
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5. 


These  may  be  looked  upon  as  "assured"  values  of  gain--the  array  will 
yield  at  least  this  much  at  this  frequency  at  any  azimuth  and  despite 
an  error  in  sound  speed  as  great  as  that  considered  (+  100  ft/sec). 
These  "assured"  values  fall  short  of  the  "asymptotic"  values  (6.02, 
9.54,  and  12.04  db)  by  an  amount  which  becomes  progressively  larger 
as  the  number  oi  sensors  is  increased.  Thus,  for  example,  quadrupling 
the  number  of  sensors  from  4  to  16  yields  only  4.38  db  instead  of 
the  6.02  db  which  might  have  been  hoped  for.  Clearly  the  gain  tends 
to  improve  slowly  with  increase  in  the  number  of  sensors;  unless 
steps  are  taken  to  measure  the  sound  speed,  and  thus  to  reduce  the 
major  cause  of  loss  of  gain,  the  maximum  practical  number  of  sensors 
to  include  in  such  an  array  probably  does  not  exceed  16. 

However,  the  negative  aspect  of  the  foregoing  remark  should  not 
obscure  the  salient  fact  that  such  arrays  offer  appreciable  gain--up 
to  almost  10  db--and  this  much  gain  implies  a  very  significant 
improvement  in  the  range  at  which  a  given  yield  can  be  detected. 
Furthermore,  the  gain  rises  with  frequency--the  direction  whtch 
favors  the  detection  of  small  yields. 
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VI.  DIRECTIVITY  PATTERNS 

The  purpose  of  employing  arrays  such  as  these  is  to  achieve 
signal- to-noise  ratio  improvement  (for  both  turbulent  and  acoustic 

•w 

noise)  so  as  to  improve  the  detection  of  weak  signals.  It  is  not 
the  purpose  to  use  such  phased  arrays  to  determine  the  azimuth  of 
the  signal;  indeed  these  arrays  would  be  relatively  poor  for  this 
purpose,  and  azimuth  usually  can  be  determined  oetter  by  correlation 
methods.  However,  the  directivity  patterns  of  these  arrays  are  of 
some  interest,  not  only  to  see  why  they  are  not  very  good  for  azimutn 
determination,  but  also  in  order  to  judge  how  many  phasing  azimuths 
are  needed  to  cover  2tt  radians. 

Figures  5-8  show  sample  directivity  patterns  of  the  continuously 
sensitive  array  and  of  the  4-,  9-,  and  16-sensor  arrays.  In  these 
figures  no  error  in  sound  speed  is  introduced;  such  an  error  would 
lead  to  even  more  irregular  patterns,  especially  at  the  higher 
frequencies.  For  each  array,  nine  patterns  are  shown,  for  phasing 
co  0°,  22.5°,  and  45°  azimuths,  and  for  frequencies  0.3,  1.0,  1.5  cps. 
Ordinarily  a  frequency  of  1.0  cps  is  regarded  as  the  upper  end  of 
the  band  of  interest;  1 .5  cps  is  included  in  these  figures  as  a 
matter  of  interest  to  show  how  the  patterns  behave  at  higher  fre¬ 
quency.  It  will  be  seen  that  side  lobe  heights  vary  with  both 
frequency  and  steering  azimuth. 

As  antenna  patterns  go,  these  are  poor;  the  side  lobes  are 
numerous  and  high--many  about  as  high  as  the  mam  beam.  The  main 
beam  is  relatively  blunt  and  not  capable  of  very  precise  azimuth 
determination,  even  apart  from  the  possible  ambiguities  of  side  lobe 
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reception.  Perhaps  it  should  he  emphasized  that,  poor  though  these 
he,  they  do  exhibit  gain. 

Tnasmuch  as  an  array  which  is  phased  to  some  particular  azimuth 
can  "cover"  adequately  only  a  modest  range  of  angles  to  each  side,  it 
is  necessary  to  arrange,  one  way  or  another,  to  phase  the  array  to  a 
variety  of  azimuths  in  order  to  cover  the  full  2n  radians.  This  night 
be  done  by  looking,  in  time  sequence,  in  one  direction  after  another, 
or  it  might  be  done  by  looking  simultaneously  in  numerous  directions. 
The  merits  and  drawbacks  of  these  methods  are  not  part  of  this  dis¬ 
cussion;  in  any  case  it  is  important  to  know  in  how  many  different 
directions  it  is  necessary  to  look. 

Figure  S  shows  the  dependence  of  the  main  beam  crossover  points 
upon  number  of  beams  used  at  a  frequency  of  1  cps  (the  customary 
upper  frequency  limit).  These  curves  are  based  upon  the  beamwidth 
for  the  array  phased  to  azimuth  ■  0°;  the  actual  crossover  points 
would  differ  slightly  from  beam  to  beam  because  of  the  slight  change 
in  beamwidth  with  azimuth.  Furthermore,  Fig.  9  assumes  no  error  in 
sound  speed;  an  error  would  alter  these  values,  perhaps  appreciably 
at  certain  azimuths.  However,  this  figure  provides  a  sufficient 
basis  upon  which  to  estimate  the  number  of  beams  needed.  The  number 
needed  depends  upon  the  crossover  value  which  is  accepted  and  upon 
the  upper  frequency  to  b*.  covered;  tabulating  the  number  of  beams 
needed  versus  the  crossover  point  at  1  cps: 
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Number 

of  Reams  to 

Cove  r 

2tt 

4  sensors 

9  sensors 

16 

sensors 

23 

18 

17 

25 

20 

18 

27 

22 

20 

30 

24 

22 

33 

27 

24 

36 

29 

27 

Clearly,  the  16-sensor  array  is  slightly  advantageous  in  this  respect 
because  its  beamwidth  is  appreciably  wider  than  the  beamwidths  of  the 
sparser  arrays.  This  advantage,  together  with  the  somewhat  better 
gain,  probably  justifies  the  added  expense  of  providing  16  sensors 
rather  than  9.  One  can  infer,  from  the  figure,  about  how  many  beams 
would  be  needed  for  arrays  of  still  more  sensors,  or  for  arrays  of 
different  size. 
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VII.  CONCLUSIONS 

The  use  of  phased  arrays  of  omnidirectional  acoustic  sensors 
(which  themselves  contain  devices  to  reduce  turbulent  noise)  would 
provide  significant  improvement  of  both  the  signal-to-turbulent- 
noise  ratio  and  the  signal-to-acoustic-noise  ratio  as  compared  with 
the  output  of  a  single  such  sensor. 

The  amount  of  improvement  is  greatest  at  the  upper  end  of  the 
frequency  band--thus  tending  to  improve  the  detection  cf  weaker 
signals  from  smaller  yields  relative  to  large-yield  signals.  A 
16-element  array  offers,  at  its  best  frequency,  about  10  db  or  more 
of  gain  at  any  azimuth  even  without  correcting  for  sound  speed  errors 
of  about  +  100  ft/sec. 

To  cover  the  full  2tt  radians  at  1  cps,  a  1/2  x  1/2  mile  16-element 
array  necessitates  looking  rt  about  22  azimuths  (for  -3  db  crossover)  . 

The  directivity  patterns  are  rather  poor  for  azimuth  determina¬ 
tion.  Thus  the  utility  of  such  an  array  would  stem  from  the  pro¬ 
nounced  improvement  in  sensitivity  for  first  detection  of  weak  signals, 
not  in  improved  azimuth  determination. 
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DIRECTIVITY  PATTERNS  AND  ARRAY  GAINS 


In  this  appendix  the  expressions  for  the  directivity  patterns 

•k 

and  array  gains  discussed  in  the  text  are  listed  without  derivation. 
The  notation  defined  in  the  text  is  used. 

Let  one  edge  of  the  square  array  be  of  length  L. 

Let 

9  =  2TTfL/c 

8  =  2rrfL/c 

o  o 

Then  the  directivity  patterns  (amplitude)  for  these  arrays  are: 
for  N  =  A 

S  =  2  |  cos  Tf  (cos  qj  +  sin  <p) 

+  cos  |j|-  (cos  cp  sin  cp) 


The  expressions  given  assume  that  the  effective  acoustic  loca¬ 
tion  of  each  sensor  does  not  shift  with  frequency.  Inasmuch  as 
each  sensor  has  a  noise-cancelling  pipe  array  attached  to  it,  this 
requires  that  the  pipes  be  symmetrically  disposed  about  each  acoustic 
sensor.  This,  in  turn,  dictates  that  the  pipes  extend  somewhat  beyond 
the  boundary  of  the  jl / 2  x  1/2  mile  region.  However,  the  relative 
insensitivity  of  the  array  gain  to  errors  in  c  suggests  thet  this 
requirement  of  symmetry  is  more  a  formality  than  a  practical  necessity. 


8  1 

-  ~  (cos  0  +  sin  0)  J 

-  ~  (cos  0  -  sin  0)  "I  | 
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for  N  -  9 


r  ro  0 

1  +  2  |  cos  |  (cos  cp  +  sm  9)  -  (cos  ft  +  sin  6)  I 


+  cos 


I 


2  (cos  9  -  sin  9)  -  ~  (cos  0  -  sin  0) 


+  COS  cos  Cp  ■  2“  cos  e 


+  cos 


S  Sir. 


<P 


sin  0 


]} 


for  N  -  16 


S  =  2  i  cos 


{  cos  |f  <cos  «P  +  3i-n  9)  -  ~  (cos  0  +  sin  0)  J 

re  0 

+  cos  II  (cos  9  -  sin  cp)  -  (cos  0  -  sin  0)  J 
_  ^ 

+  cos  [*■  (cos  9  +  j  sin  9)  -  (cos  9  +  «  sin  0)  J 

+  cos  [f  (cos  9  -  j  Sin  9)  -  (cos  0  -  j  sin  6)  J 

P 

+  cos  [f  (3  cos  <P  +  sin  9)  -  ~~  (1  cos  0  +  sin  0)  1 


+  cos 


ffi  1  3  1 

||  (-  cos  9  -  Sin  9)  -  ~  (|  cos  0  -  sin  0) 


fB  .1  1 

+  cos  ||  ("3  cos  cr  +  ~  si>  9) 

+  cos  P|  (~  cos  9  -  ~  sin  9) 


p  1  1 

~  (-j  cos  0  +  ~  sin  9)  ! 

8o  1  1 

(-  cos  0  -  -  sin  0)  "j  j 
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for  N  =  * 


S  =< 


0, 


5 in  cos  qj  _  _£  cos  q  } 

[f  cos  cp  -  —  cos  e] 


As  to  array  gain,  the  numerator  of  the  gain  expression  poses  no 
problem;  it  is  only  necessary  to  set  tp  =  0  in  the  foregoing  pattern 
expressions,  and  then  precede  to  evaluate  the  numerator.  The  difficulty 
resides  in  the  evaluation  of  the  denominator.  Let  the  denominator  be 

2tt 

D-^j  S2d, 
o 

Then,  for  these  arrays: 
for  N  =  4 

D  =  4  +  4  J  (3)  f  cos  [3  cos  Q]  +  cos  [9  sin  9]} 
o  o  o 

+  2  Jq(/2  9)  {  cos  [/2  9o  cos  -  9)]  +  cos  [fl  9q  sin  -  0)]} 

for  N  ■  9 

D  *  9  +  6  J^(3)  |  cos  f3o  cos  0I  +  cos  j^90  sin  9*1  j- 

+  12  Jq  (  ^  )  i  COS  jjp  cos  9*j  +  cos  r~  sin  0*]  } 

+  2  Jc  ^/2  9^  j  cos  |”/2  3o  cos  -  9)1  +  cos  f/2  0o  sin  -  9)  1 

+  8  Jo  (2^  p)  {  cos  I  2^  Bo  COS  (4  "  e)1'  +  COS  Bo  Sin  (4  "  6>  1 


(continued) 
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+  4  J 


(*F  P)  {  COS  1.4  Po  COS  (crl  '  e>  1  +  cos  [4  30  (, 0fl  - 


9) 


+  cos 


[V*  Bo  CO"  ^®l  +  0)]  +  cos  [4  3o  sin  (flfj  +  9) 


where  sin  a 


1  'A 


cos  ai  =  A 


for  N  *  16 


D  =  16  +  8  Jo  (g)  {  cos  cos  9*|  +  cos  sin  9 

+  ^  ^o  (  3  B)  {  COS  I  3  So  COS  +  cos  £3  *0  sin  sj  j- 

+  24  Jo  (  3  S)  {  COS  Li  eo  COS  e]  +  cos  [i  S„  sin  e]  } 

+  2  Jo  (//2  13 )  {  cos  [/*  B0  cos  (f  -  *>)]  +  cos  [,/T  9o  sin  (|  -  9)^  j 

+  8  Jo  (¥  8)  {  ‘°S  Bo  cos  (|  -  9)]  *  cos  [*£  B0  Sin  (2  -  9)1  } 

+  18  Jo  (  f  B)  {  “>s  [f  f0  cos  (f  -  9)1  +  cos  [2Z  Bo  sin  ("  9)1  1 


+  6  J 


O  (4^  B)  {  C0S  [4^  Po  COS  (ar2  "  e)]  +  COS  [4F  Po  sin 


-  0) 


+  cos  [4^  3o  cos  (a2  +  9)]  +  cos  [4^  eo  Sin  (c*2  +  9)1  } 
12  Jo  (4^B)  {  COS  [4  So  COS  (ai  "  6>]  +  cos  r4 


+  cos  y-  Po  sin  (or  -  9)  I 


r./5  a 

+  cos  I  3  cos 
i  3  o 


(»!  +  9)J  +  cos  M  Bc  sin  +  9)  |  J 

4  Jo  (  4^  S)  {  cos  [41  Bo  C0S  (a3  -  £)  :  +  e<>s  |4^  Bo  Sin  (o3  -  9)  1 

P  /l3  r/f^ 

+  cos  }«-—  3o  cos  ( Of 3  +  9)  +  cos  I  *4  9o  sin  (a^  +  9)  |  \ 
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where  sin  o>2  -  cos  a2  -  7^ 

Sln  *3  =  jh  C°S  *3  =  7H 

for  N  «  °o 

No  expression  for  D  has  been  obtained.  The  points  needed  to 
plot  Fig.  1  were  obtained  by  numerical  integration  on  a  digital 
computer.  This  is  a  relatively  slow  and  expensive  routine,  and  it 
was  not  thought  to  be  worthwhile  to  explore  the  full  domain  of  c  and 
8,  as  was  done  for  the  finite  values  of  N. 


Note  that  these  expressior s  for  D  have  a  common  form 

D  =  N  +  a  sum  of  terms  of  the  form 

n  J  (kg)  {  cos  cos  (cr  -  9)] 

o  o 

+  cos  [k@o  sin  (0  -  0)]  } 

It  ’■■■•11  be  seen  that  n,  k,  and  Of  are  related  to  the  geometry  of  the 
array  itself  as  follows: 

0  Connect  every  sensor  point  in  the  array  to  every  other 
point.  In  doi'ig  so,  lines  which  are  themselves  component 
parts  of  longer  lines  are  nevertheless  to  be  regarded  as 
distinct . 

0  Then  n  ■  the  total  number  of  such  lines  whose  length  is 
kL  and  which  lie  at  azimuth  a  in  the  array. 
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For  example,  in  the  array  of  4  sensors,  there  are  tour  lines 
(the  edges)  of  length  such  that  k  =  1  and  whose  azimuths  are  either 
0  or  n/2  (and  these  combine  because  cos(tt/2  -  a)  ~  sin  a),  and  there 
are  two  lines  (the  diagonals)  whose  length  is  such  that  k  =  fT  and 
whose  azimuths  are  +rr/4  (which  also  combine).  In  the  9-element 
array  these  same  6  lines  are  seen,  plus  two  more  lines  of  length 
corresponding  to  k  =  1  passing  through  the  center,  plus  lines  with 
k  *  /l! 2,  k  =  1/2.  i  and  k  =  /T/2. 

This  structure  for  the  expression  for  D  holds  for  all  finite  N. 
The  proof  of  this  rule  is  not  difficult. 
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